Localization of DNA damage by current exchanging repair enzymes: 
effects of cooperativity on detection time 
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How DNA repair enzymes find the relatively rare sites of damage is not known in great detail. 
Recent experiments and molecular data suggest that the individual repair enzymes do not work 
independently of each other, but rather interact with each other through currents exchanged along 
DNA. A damaged site in DNA hinders this exchange and this makes it possible to quickly free up 
resources from error free stretches of DNA. Here the size of the speedup gained from this current 
exchange mechanism is calculated and the characteristic length and time scales are identified. In 
particular for Escherichia coli we estimate the speedup to be 50000/iV, where N is the number of 
repair enzymes participating in the current exchange mechanism. Even though N is not exactly 
known a speedup of order 10 is not entirely unreasonable. Furthermore upon over expression of 
repair enzymes the detection time only varies as iV -1 / 2 and not as 1/N. This behavior is of interest 
in assessing the impact of stress full and radioactive environments on individual cell mutation rates. 
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PACS numbers: 

Functional and properly working DNA is a prerequi- 
site for the survival of all cells and organisms. However 
all the time localized errors in the DNA string are intro- 
duced both by external stress factors and normal cellular 
functions like metabolism. It is thus of utmost impor- 
tance to quickly detect the damage and repair the DNA. 
In Escherichia Coli an enzyme MutY is responsible for 
localizing 'oxidative damage' to single base pairs. Ho- 
mologous enzymes presumably with a similar mode of 
action is found in most other species. It is well known 
that MutY is able to slowly move along the DNA Q 
while it checks the integrity of each base pair. However it 
has previously been pointed out that MutY acts surpris- 
ingly fast 0]. In fact it probably localizes a defect base 
pair faster than is possible by the slow scanning above. 
Based on microscopic data it has recently been hypothe- 
sized that two MutY complexes bound to DNA are able 
to communicate via currents in the DNA and in this way 
speedup the localization process ■ More specifically an 
error free stretch of DNA is a good conductor, while a 
defect base pair introduces a huge resistance |4|; thus if 
a MutY enzyme receives an electron from an upstream 
MutY enzyme it knows the stretch of DNA ahead of it 
is error free and it then detach from the DNA and in- 
stead find another stretch of DNA to scan. Intuitively 
it is clear this quick freeing up of MutY enzymes from 
error free stretches of DNA speeds up the localization of 
damaged base pairs, but by how much? There are two 
relevant time scales in the proposed process. The first is 
the time T it on average takes to localize a damaged base 
pair by slowly scanning the DNA, without utilizing the 
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currents. The second is the time r it takes to realize a 
stretch of DNA is without errors using the current, i.e. r 
is the time from a MutY enzyme attach to an error free 
piece of DNA until it detaches and goes somewhere else to 
attach. In this paper, we first show that the time it takes 
for MutY in the above scheme to localize a damaged base 
p air is roughly \fTr corresponding to a speedup of order 
y/T /t. This expression for the speedup is even retained 
in the presence of many other kinds of current interacting 
repair enzymes. However T is in this case the time until 
the first repair enzyme, which need not be MutY, locates 
the error by scanning alone. 



I. MODEL 

The model is pictorially presented in Figure ^ See 
also The repair enzyme MutY contains an evolu- 

tionary well-conserved [4Fe4S] cluster that is suspected 
to change its charge configuration from +2 to +3 upon 
binding to DNA .3}. Upon binding to DNA an electron 
is thus emitted into DNA, while upon receipt of an elec- 
tron from DNA the MutY-DNA binding is destabilized. 
As only error-free stretches of DNA are able to trans- 
port the electron from one MutY enzyme to a neigh- 
boring MutY enzyme this charge exchange enables 
MuTY to quickly free up scanning resources from error 
free stretches of DNA Q — see Figure In the original 
proposal only MutY enzymes participates in the current 
exchange and to illustrate the line of thought we first 
consider this scenario. Afterward, we then extend the 
consideration to the case where there are many different 
kinds of repair enzymes each specialized to fix a specific 
kind of damage. Furthermore, we also consider the effect 
of a finite scan length before MutY spontaneously drops 
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FIG. 1: a) The left MutY repair enzyme is bound to DNA 
and slowly progress to the right while it scans the integrity 
of the base pairing. The [4Fe4S] cluster in MutY is in a +3 
charge configuration when bound to DNA, while it is in a +2 
configuration when not bound (right MutY). b) Upon binding 
to DNA the right MutY enzyme emit an electron onto DNA 
and changes the charge of its [4Fe4S] cluster from +2 to +3. c) 
If the DNA is free of errors the emitted electron travels along 
the DNA until it reaches the left MutY enzyme. Here the 
electron changes the charge of the [4Fe4S] cluster to +2 and 
thus destabilizes the DNA binding of the left MutY enzyme. 
The left MutY enzyme then attach to and scan a different 
section of DNA more likely to posses an error, d) If on the 
other hand the DNA piece in between the two MutY enzymes 
contains an error the electron never reaches the left MutY 
enzyme, which then keep scanning the DNA until it reaches 
and fixes the error. The current exchange thus selectively free 
up resources from error free patches of DNA. 



off DNA. 



A. Only MutY and it scans forever 

After a a base pair get damaged, how long time does 
it take before a MutY enzyme has located the error? In 
the regime where current scanning yields a considerable 
speedup in the location of the error, the damaged base 
pair typically is located by a MutY enzyme that binds to 
DNA next to it and then scans all the way up to and in- 
cluding the faulty base pair (and not by a MutY enzyme 
that just happened to be nearby at the time of damage). 
Let tdetect denote the typical detection time and v the 
scan velocity of MutY. The MutY enzyme then typically 
docks onto DNA within a distance of the order ^detect- 
The rate with which a MutY enzyme randomly docks 
onto a specific base pair and start scanning is denoted by 
k. The probability for a MutY enzyme to land within a 
distance of videtcct of the error in the time interval idotoct 



can be estimated as kvt\ etect . As at time idetect a MutY 
enzyme typically arrives at the faulty base pair this prob- 
ability is of order 1 and 
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A more detailed model calculation in the supplementary 
material yields the same result apart from a factor 1.3. 
The average docking rate k can be expressed as 



k = 



(MutY) 

71 



vTt 



(2) 



where r is the time between two successive binding events 
for a single MutY enzyme. A^ (MutY) is the total number 
of MutY enzymes and L is the total number of base pairs 
in DNA. T = L/v/N(MutY) is the time it takes for all 
the MutY enzymes to scan all the bases of DNA once. 
From standard Poisson statistics T is also the detection 
time of a damaged point in the traditional scenario with- 
out any cooperation between MutY enzymes. In terms 
of T and r the detection time is 
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The cooperation thus give a speed up of order \JtJT. 
Before I get down to consider numbers and orders of mag- 
nitude I want to make the model a tad more realistic. 



B. Full model 

The functionally central [4Fe4S] cluster is also present 
in other repair enzymes e.g. endonuclease III. It is thus 
very likely that other repair enzymes also are able to 
inject currents into DNA and participate in the electri- 
cal scanning of DNA. Furthermore all of these repair en- 
zymes then get 'attracted' to the damaged DNA pair in 
exactly the same way as MutY. In the above model and 
calculation we thus have to replace 'MutY' by 'any repair 
enzyme participating in the DNA mediated charge trans- 
port' — for short just a repair enzyme below. Likewise 
the calculated detection time ^detect becomes the time be- 
fore the first repair enzyme finds the damaged site and 
T the average time it take for any repair enzyme to find 
the site without using currents. I have here implicitly as- 
sumed that both the scan velocity v and the single repair 
enzyme attempt frequency 1 /r are of the same order of 
magnitude for all repair enzymes i.e. MutY is a typical re- 
pair enzyme. Biologically the time detect is not the most 
relevant one as the first repair enzyme that arrives at the 
damaged base pair very likely is unable to fix the damage. 
For instance is MutY's primary target oxidation damage 
resulting in G:A and 8-oxo-G:A mismatched pairs. In- 
stead the first MutY enzyme arrives on average as the 
N/N {MutY) repair enzyme resulting in a detection time 
of order A^/A^(MutY)idctoct- Here N is the total number 
of repair enzymes. N/N(MutY) can also be expressed 
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as T/T (MutY), where T(MutY) is the time it takes for 
a MutY enzyme to find the site by scanning alone. The 
MutY detection time is thus 

i dotcct (MutY) « T(MutY) , (4) 

again corresponding to a speedup of size y/rJT, but this 
time T is the time it take for any repair enzyme to locate 
the damage. 

I have not yet considered that MutY is known sponta- 
neously to drop off DNA after having scanned of the order 
100 base pairs (bp) In order to estimate the effect 
of this I am going to derive the above result Eq. Q in a 
slightly different manner. The MutY enzyme that even- 
tually localizes the damage typically docks onto DNA 
within a distance A from the faulty base pair, where A 
both has to be small enough to allow scanning all the 
way up to the error i.e. < lOObp and also so small that 
it is unlikely another repair enzyme dock in front of the 
MutY enzyme while it scans. The last le ngth scale can 
be estimated to be of the order A = \pvjk — vVTt. 
The MutY detection time idetect(MutY) is then deter- 
mined as above by setting the probability for a MutY 
enzyme to dock within a distance A in the time interval 
idetect(MutY) equal to 1 i.e. 

VT 

idctcct(MutY) = T(MutY) — , (5) 

where A = min(100 bp, vy/rT). Notice that the choice 
A = vVtT leads to Eq. 0. 

C. Estimating order of magnitude 

I have not been able to find any direct experimen- 
tal measurements of r and T, so I am instead going 
to use a somewhat more uncertain way to estimate the 
size of the reduction —. The distance vt is the aver- 

VT 

age scan length of a repair enzyme (MutY) in the pres- 
ence of the current interactions. The numerator A is 
the smallest of the maxim al scan length lOObp and the 
docking distance \JvtvT. In fact, I believe these two 
distances are of the same order of magnitude as any- 
thing else seems inefficient, so I assume VvtvT < lOObp, 
with equality as the most likely option. The distance 
vT = L/N is the average distance between repair en- 
zymes. From this I can estimate vt. The reduction is 
thus A > vT/100 bp = 5 • 10 4 /Y, where N is the total 
number of repair enzymes with a current exchange mech- 
anism similar to MutY and I have used that the length 
of E.Coli's DNA is 5 • 10 6 base pairs. Unfortunately N 
is unknown. In |{| the number of the two [4Fe4S] 2+ con- 
taining repair enzymes MutY and endonuclease III is es- 
timated to be 30 and 400 respectively. In the same paper 



the number of formamidopyrimidine glycosylase (FAPy 
or MutM) repair enzymes is estimated to 400. FAPy does 
not contain the [4Fe4S] 2+ cluster. The target of FAPy is 
8oxoG which is estimated to constitute 5% of all adducts 
due to oxidative damage . All in all it seems reasonable 
that the total number of repair enzymes participating in 
the current exchange mechanism is significantly smaller 
than 50000, and that a speedup of order 10 is realistic. 
Notice this would correspond to a typical scan length 10 
times smaller than the maximal one (100 bp) due to the 
current exchange mechanism (vt is of the order 10 bp). 



II. CONCLUSION 

In this paper we have considered some of the implica- 
tions of a very interesting proposal for cooperation be- 
tween repair enzymes in the localization of defects in sin- 
gle base pairs. First we have pointed out that the mech- 
anism is likely to speedup the localization by a factor 
of order 10 compared with independent scanning by the 
repair enzymes. If an error is detected by the 30 MutY 
enzymes in say 20 minutes (it has to be considerably 
shorter than the replication time 0) this corresponds 
to a reduction in scanning speed from 125 bp/sec to 13 
bp/sec. For comparison, the scan velocity for RNA poly- 
merase is 50 bp/sec, while for DNA polymerase it is 1000 
bp/sec. Another in principle testable prediction is that 
upon over expression of all the repair enzymes with say a 
factor 4 the detection time is only decreased by a factor 
2 (both T(MutY) and T are 4 times smaller), because 
information is not carried as efficiently along the DNA 
any more. Physiologically oxidative and radiative envi- 
ronments may result in an increased expression of repair 
enzymes 0, so the square root behavior and the cou- 
pling of the effectiveness of different kinds of repair en- 
zymes is potentially of huge importance for the mutation 
rates in these kinds of stress full environments. Summing 
up, only further experimentation can finally confirm this 
charge transport mechanism, while we here have demon- 
strated that the mechanism indeed offers a great benefit 
for the cell. In addition, the model is a nice toy model 
for protein cooperativity and one might wonder if the 
underlying principles behind could be of practical use in 
apparently unrelated engineering problems. 
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